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CHAPTER I 


INTRODUCTION 


1 . 1 Origin and Importance of this Investigation 

Cavitation is an important consideration in the design and 

^ i ^ ^ ■ ^ : . " 

analysis of liquid handling machinery. Cavitation will occur in any 

liquid flow whenever the local pressure is reduced below some critical 

'i - ' ' " 1 ' 

value It is generally assumed that cavitation will occur when the 
minimum pressure in a flow reaches the vapor pressure based on the 
free scream temperature of the liquid. However, this assumption is 
invalid if significant thermodynamic effects are present. 

A continuous supply of vapor is required to sustain a deveicped. 
cavlcy Ihis vaporization process requires a transfer of heat to the 
cavity, so that the temperature of the cavity will always be less , 
chan the free st'eam temperature of the liquid. This localized 
cooling phenomenon is referred to as the thermodynamic effect. 

The determination of the cavity pressure is of primary 
Importance in developed cavity flows. The thermodynamic-ef feet can 
significantly influence the cavity pressure when the liquid 
temperature (T) is close to the critical temperature where the 
magnitude of the vapor pressure (P^) and its derivative dP^/dT are 
large. i , 

/. 
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In recenc years, thermodynamic effects have' become Increasingly 
Important due to the need to design turbomachines for operation in 
cryogens such as liquid hydrogen and liquid oxygen. For these liquids, 
the operating temperature is closer to the critical temperature, and 
thermodynamic effects can be signlif leant, 

A significant thermodynamic effect can occur in any cavitating 
liquid Due to the complex nature of a developed cavity, prediction 
of cavity temperature depressions is a difficult task, and extensive 
study of this subject is needed^ 

1 _ 2 Previous Investigations i 

5^ ■ i j ; : : ' j 

Stahl and Stepanoff [1] were two of the first Investigarors 

i. : i : 

who studied the thermodynamic effect on developed cavities : In that 

paper, the B-factor method based on a quasi-static theory was 

I ; 

developed tor predicting thermodynamic effects with particular 

I I : 

emphasis on pump applications. Fishhr [2], Jacobs [3], and Acosta 
and Hollander [4J also considered the B-faetor approach and its 
applicat ion Later , other Investigators such as Salemann [5] , i 
Spraker [6], and Hammict [7] experimentally determined the effect 
of fluid property variations on the performance characteristics of 
pumps - ' i ' I i ! 

NASA has undertaken a program to extend the previous work in 
an effort to obtain improved design criteria to aid in the 

Numbers in brackets refer to the List of References at the end of 
this thesis. 
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predlct-ion of cavitating pump performanceo Ruggeri, eto al. [8-10] 
and Hord, et al. [11-16] have obtained; extensive experimental data 
of cavity pressure and temperature depressions for a variety of 
model shapes and fluids, and have correlated the results using the 

B-factor theory , ; j 

Recently, Billet [17] has proposed an alternative to the 
B-factor approach This method, called the ehtrainmentr^cheoryi Is j 
a dynamic approach based on an energy balance tor the cavity ^ The 
entrainment theory appears to be a better physical model than the 
B-factor theory, although both theories are quasi-empirical in 
nature 

1 5 Purpose cf th:s Investigation 

The purpose of this investigation is twofold, Flrsjt, the 
entrainment theory concepts are extended by additional experimehtal 
and theoretical anaiyses to show chat It; is a favor able alternative 
tc the B-taccot theory. Second, a detailed summary or other 
investigations is given tc establish the current state of knowledge 
for predicting the ,-hetmodynamic effect on developed cavitation. 

Ihe principal tbe?r.sticai and experimental results are presented 
fcr comparison with results obtained in this study 


CHAPTER II 


EXPERIMENTAL INVESTIGATIONS 

2 . 1 Introduction 

Due to the complex nature of a developed cav/icy flow, experi- 
mental determination of the cavity characteristics is probably the 
most important phase of this investigation. Knowledge of the cavity 
pressure, temperature, geometry, and vapor mass flow rate is 
necessary to treat problems concerning thermodynamic effects as 
well as other developed cavity phenomenon. In this chapter, 
experimental procedures used in this study are presentecl and the 
chara: teristics of the resulting data discussed 

2 . 2 Facilities and Models 

Two facilities were used to conduct the experimental 
invest igation, both of which are located at the Applied Research 
Laboratory of The Pennsylvania State Universityc. The primary 
faciiity used is the NASA-sponsored uicrn-high-speed cavitation 
tunnel as shewn in Figure 1-. The tunnel is a recirculating system 
having a. circular 3.8 cm internal diameter test section, and has 
the capability of operating over a viide range of velocities, 
pressures, and temperatures in various liquids The working fluids 
used in this study are water and Freon 113. A detailed description 

H- 
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of the facility and data acquisition system jisi given in References 
[17-22]- A second facility, a 30.5 cm diameter test section water 
tunnel with a more limited operating range was used for the 
venrllated cavity tests,- A detailed description of this facility 
is given in References [22-23], 

A total of fourteen sting-mounted ogive test models having two 
basic nose contours were employed in the exper-imeiits : Thte zero- 

== ' I : ! 

caliber ogives have a blunt nose, whereas the q^uarter-caliber ogives 
have a rounded nose with a radius of curvature equal to one quarter 
of the model diameter, A description of the two basic nose contours 
is given in Figure 2, and a summary of the test models used is given 
in Table I - , , 

The six mcdels that were used for the ven,t.ilated cavity tests 
were three zero-caliber ogives with model diameters of 0-318, 0-635, 
and 1 27 cm, and three q^uaijter-caliber ggives with model diameters i 
of 0 318, 0635, and 1.27 cm. These models have a hollow center 
from which air is injected through holes near the leading edge ho 
form the ventilated cavities, and a tube along the surface of the 
model with a pressure port close to the leading edge to measure the 
ca’' ify pressure- 

Four models were used to measure the cavity pressure of natural 
developed cavitleis:i two zero-caliber ogives having^-mpdel diameters 
of 0,318 and 0- 635 cm, andfitwo quartef^'caliber ogivesc having model 
diameters cf 0,318 and 0 6.35 cm. These modaii harve three tubes 
along the surface of the model with each tube having a pressure 



L...J 

pert located at a different axibl position on the model, allowing the 
axial variation in cavity pressure to be determined, A photograph 
of the four cavity pressur'e_ iiodels is given in Figure 3. 

Four models were used to measure the cavity temperature for 

t I 

natural developed cavities as shown in Fi gure 4 The two' zero-caliber 
ogives have model diameters of 0.318 and 0,635 cm and the two quarter- 
caliber ogives also have model diameters of 0,318 and 0»635 cm. 

These models have three ports in which thermocouple beads are 
cemented, and the thermocouple leads exit the tunnel tlirough the 
hollow center of the model an,d the sting mount. The thermocouples 
are mounted at three different axial positions on the model so that 
the axial distribution of temperature within the cavity could be ' 
determined Details of the thermocouple fabrication i and installation 
procedure are given in Appendix A. In addition, the two larger 
diameter models have one tube along the surface of the model to 
mcnitor the cavity pressure. 


2 , 3 Cavity Geometry , . 

The geometric characteristics of developed cavities were 


profiles taken in - 


determined from photographs of developed cavity 
both the 30 5 cm and 3,8 cm tunnels for both model contours. 

Typical photographs of developed cavities for the two model contours 
are shown in Figures 5 and 6 for water and Figures 7 and 8 for 
Freon 113, From measurements of the cavity profile taken from the 


photographs, the maximum cavity diameter (Dj^) » cavity half length 
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(A) , and cavity surface area (A^) were determined for both natural 
and ventilated cavities at various velocities, dimensionless cavity 
lengths, and model diameters. Details of the procedure and resulting 
data are given by Billet [17] and Billet, Holl, and Weir [24] 

The data show that for a given model geometry and cavitation 

! ; 

number, ventiiared and natural cavities have the same profile, within 
experimental error, independent of the free stream velocity, and that 
the ravity dimensions are directly proportional to the model diameter 
This is in agreement with Waid [25], who measured developed cavity 
profiles behind a circular disk, In addltlcn, for the same flow 
conditions, the surface area of a cavity on a zero-caliber ogive 
IS about twice that of a quarter-caliber cgive The resulting data 
ter rhe area coefficient (A^/D^) as a function cf dimensionless 
C3'’i y length (I'D) for the two model contours are given in Figures 
9 and lO 

2 4 Cavity Flew Coefficients 

The cavity flow coefficient (C^) is a dimensionless representa- 
tion cf the volume flow rate of gas required to sustain a ventilated 
cavity and is defined as 


In general, the flow coefficient is a function of model geometry, 
free stream velocity, and cavitation number. 


= _ 2 _ 


V D" 


( 1 ) 


Measurements were made in the 30,5 cm water tunnel of the volume 


flow rate of gas and the cavitation number at various velocities and 
dimensionless cavity lengths for the six ventilated cavity models 
The gas flow rate was measured using a Gilmont float-type flowmeter 
The free stream pressure was adjusted so that the cav/ity pressure was 
equal to the partial pressure of gas in the free stream to ensure that 
minimal gas diffusion occurred across the cavity wall From the data, 
the flow coefficient (C^) as a function of cavitation number (a) , 
and the cavitation number as a function of dimensionless cavity 
length CL/D) were determined for various flow conditions. Details of 
the experimental procedure are given by Billet and Weir [26-27], 

The resulting data of the flow coefficient as a function of 
cavitation number is presented in Figures 11 to 13 for the three 
zero-caliber ogive test models and in Figures 14 to 16 for the three 
quarter -caliber ogive test models. In addition. Figures 17 and 18 
give the cavitation number as a function of dimensionless cavity 
length for the two model contours . 

The results demonstrate that for a given model geometry, the 
flow coerticient increases with increasing velocity and model 
dlamexer and Increases with decreasing cavitation number The 
catfitacion number for a given model contour is a function of 
dimensionless cavity length only- Finally, for the same flow 
conditions, the flow coefficient for the zero-caliber ogive is 
about five times that for the quarter-caliber ogive. 
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2 o5 Cavity Pressure Measurements 

Measurements of cavity pressures were made in the 3.8 cm tunnel 
for natural developed cavities In Freon 113 and water. Pressures 
were measured using Pace variable reluctance transducers with an 
accuracy of ± 0,^07 N/cm^ and an integrating digital voltmeter to 
obtain time averaged data. The cavity lengths were determined by 
aligning the cavity closure point with lines scribed on the models. 
Cavitation numbers based on local cavity pressure at three different 
axial positions within the cavity were obtained in Freon 113 for the 
four cavity pressure models over a velocity range of 19.5 to 

36.6 m/sec., temperature range of 20 to 93°C, and dimensionless 

i 

cavity length range of 2.0 to 7.0. In addition, measurements of 
cavitation numbers were made in water during the cavity temperature 
tests with the single cavity pressure port on the two 0.635 cm 
diameter cavity temperature models^ for a velocity range of 19.5 to 
36.6 m/sec., temperature range of 60 to 120“C, and dimensionless 
cavity length range of 2.0 to 5.0. 

To determine the variation of cavitation number with 
temperature, the cavitation number measured at the pressure port 
located closest ' to the leading edge of the model was determined as a 
function of free stream temperature, velocity, and dimensionless 
cavity length. Results for the zero-caliber ogives are given in 

Figures 19 and 20 and results for the quarter-caliber ogives are 

i, 

given in Figures 21 and 22. Each data point shown is an average 
of at least five test runs. Data obtained by Billet [17] and 
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Buddenbaum [28] for Che 0o635 cm zero-caliber ogive are presented 
for comparison. Referring to the figures, ir can be seen that the 
cavitation number is independent of velocity and free stream 
temperature within experimental error, and depends oiily on 
dim(;nsionless cavity length for a given model contour and diameter. 
The cavitatioir numbers are referred to as uncorrected since the 
effect of blockage has not been tak.en into account c This effect i 
will be discussed in detail in Chapter IV, 

Cavitation numbers measured at the other axial pressure ports 
were also independent of free stream temperature and velocity. 
Figures 23 to 26 show the cavitation number as a function of 
dimeiisionless axial position (x/L) and dimensionless cavity length 
in Frebh' 113 for the two model contours. Each data point is an 
average over all temperatures and represents at least 40 test runs. 
It can be seen from the figures that the cavitation number, and the 
corresponding cavity pressure, remains roughly constant over the 
first 6Q percent of the cavity, then the cavity pressure begins to 

I , 

increase The measurements were made with the Freon 113 nearly ^ ' 
saturated with air, therefore the data is affected by both 
thermodynamic and diffused air effects. Since these effects 
cannot be separated, direct comparisons of measured cavity pressures 
and vapor pressures based on measured cavity temperatures were not 


made 
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2 6 Cavity lemperature Measurements 

Measurements of cavity temperatures were made in the 3t,8 cm 
tunnel for natural developed cavities in Freon 113 and water with 

the four cavity temperature models. Cavity temperatures were I 

1 

measured at three different axial positions within the cavity with 
0 010 cm diameter copper-constantan thermocouples with an accuracy 
of X 0=28°C, The cavity thermocouples were each connected in series 
with a downstream thermocouple so that the temperature depression (AT) 
could be measured directly = The free stream temperature was 
measured Independently with a thermocouple referenced to a O'C ice 
bath » 

All temperature readings weire taken with an integrating digital 
voltmeter to time average any temperature fluctuations c This differs 
from the procedure of Billet [17], who used a galvanometer to take 
instantaneous readings, and then only considered the minimum 
measured cavity temperatures. The time averaging technique therefore 
produces temperature depressions which are smaller than those 
obtained by Billet: but which are. more consistent with the steady- 
state entrainment analysis. 

In order to minimize the effects of variations in the amdunt 
of ncnrondensable gas dissolved in the liquids, all cavity 
temperature tests were run at nominal air contents of 14 ppm by 
mcles in water and 1200 ppm by moles in Freon 113 i However, it 
has been shown by Fricks [18] that the gas content of the liquid 
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has litcle effect on measured cavity temperature depressions for 
the models and test conditions employed in this investigation. 

Figures 27 to 29 show typical results of the measured 

temperature depression as a function of dimensionless axial 

position (x/L) r. Data obtained by Fricks [18] in Freon 113 for the 

two 0.318 cm models are presented for comparison As can be seen 

from the figures, the axial variation of the temperature depression 

along the cavity was found to be roughly linear , with the maximum 

temperature depression occurring near the leading edge of the 

cavity This is in agreement with other investigators [8-11, 13-18] 

Therefore, to consistently determine the maximum temperature 

depression (AT = T - Tp . ), the axial temperature depression 

distribution was extrapolated to the cavity leading edge as shown 

in Figures 27 to 29, and the resulting AT was considered in all 

max 

cases to be the maximum temperature depression- 

The maximum temperature depression was determined as a function 


of free stream temperature for various flow conditions Data were 
obtained in Freon 113 ever a velocity range of 19-5 to 36=6 m/sec., 
temperature range of 35 to 95°C, and dimensionless cavity length 
range of 2.0 to 7 0= In addition, data were obtained in water for 
the two 0=635 cm diameter models over a velocity range of 19.5 to 
36 6 m/sec., temperature range of 60 to l25'C, and dimensipnless 
cavity length range of 2^0 to 5.0= The results are given in 
Figures 30 to 35, where each data point is the average of at least 


10 test runs 
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It can be seen, from the figures that the maximum temperature 
depressions are about an order of magnitude larger in Freon 113 than 
in water over the temperature ranges tested. This can be attributed 

■ i ; I ' ' ' ' ^ 

to the fact that since the critical temperature of Freon 113 is 
214. I'C compared to 374, 1°C for water, the Freon 113 data were taken 
at temperatures closer to its_ critical temperature. The largest 

I 

maximum temperature depression observed was 9oO'C for the 0.635 cm 
quarter-caliber ogive in Freon 113. 

As demonstrated by Figures 30 to 35, the maximum temperature 
depression increases with increasing free stream temperature and 
dimensionless cavity length for both model contours.; However,, the 
maximum temperature depression increases with increasing velocity 
and model diameter for the quarter-caliber ogives, but decreases 
with increasing velocity and model diameter for the zero-caliber 
ogives. These results suggest fundamental differences in the 
thermodynamic characteristics for the two model contours, but this 
will be clarified by examining the data within the context of the 
entrainment theory in Chapter IV. 


2 1 Ocher Experimental Investigations 

Measurements of cavity pressures and temperatures have been 
made by other Investigators for a variety of model geometries and 
liquids These investigations are summarized in Table II, which 
gives the models and fluids tested as well as the range of flow 
parameters. The measured AT at the leading edge thermocouple is 
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considered in all cases. The data for the AT obtained in this 

ma x 

study are presented for comparison. 

The geometries of the zero- and quarter-caliber ogives have 
already been discussed. The venturi models are axisymmetrlc and 
have a quarter-round ridge at the minimum area point to ensure 
that a well defined cavity leading edge i is formed , The hydrofoil 
model is two-dimensional with a leading edge having a semi-circular 
cross-section and a streamlined trailing edge section- For a 
detailed description of the models, refer to the references cited 


in Table II. 


CHAPTER III 


THEORETICAL CONSIDERATIONS 

3 I Intr oducLioa 

A direct theoretical description of the thermodynamic processes 

involved in the maintenance of a natural developed cavity Is an 

! 

extremely complex problem It would include a detailed knowledge 
:i ’■he '.aporizatlon, condensation, convection, and diffusion 
pfo'.ejses cccuiring. Since the flow field is' two -phase, 
nonisentropic, and in most cases, turbulent, an analytical 
determination of the velocity and temperature profiles within the 
cavity and suitcundlng liquid, the ’heat required to term the cavity ^ 
dod the magnitude of any thermodynamic effects is beyond the 
current state of knowledge. 

Several simplified theories: have been developed to model the; j 
chetrncdyriamic aspects of a natural developed cavity, and are presented 
1 C the following sections. First, some preliminary considerations 
will be discussed, 

A developed cavity is described by the cavitation number (o) 
derined as 


l/2p^V^ 


0 


( 2 ) 
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in general, the cavity pressure (P ) is given by 

u 

' ‘’v + ’’g - ''■’v > «) 

where 

Py = vapor pressure based on the free stream temperature 
of the liquid (T^) , 

Pg = partial pressure of noncondensable gas in the cavity, 

and 

AP^ = vapor pressure depression due to the thermodynamic 
effect. 

Subsritution of Equation (3) into Equation (2) expresses the 
cavitation number as 


0 = 


P - 

00 


1/2P.V, 


L °o 




AP, 


1/2P V. 


L “ 


( 4 ) 


Several investigators [29-31] have shown that the contribution of 
noncondensable gas in Equation (4) can be significant, and as the 
bulk temperature of the liquid increases toward the critical 
temperature, the vapor pressure depression term will also become 
important- In addition, the vapor pressure depression can, in J 
general, be affected by the gas content of the cavity. 

The vapor pressure depression can be expressed in terms of 
the difference between the free stream and cavity temperatures. 


17 


The vapor pressure depression can be expanded in a Taylor series as 



dT 


d'^P, 


At + 


dT' 


AT' 



(5) 


where AT = T^ - T^ is the temperature depression. Now, retaining 
only the first term of the series and using the Clausius-Clapeyron 
relationship for yields 



Py A AT 

- 


1 

L- 


( 6 ) 


One can estimate the vapor pressure depression by Equation (6) once 
the temperature depression is known, provided the vapor pressure 
depression is small. 


32 Entrainment Theory 

The entrainment theory Is a semi-empirical method ±or predicting 
thermodynamic effects on developed natural cavities -i The theory, 
based on a simple energy Jalance, correlates the temperature 

. i 1 

depression (AT) in terms of the relevant nondimenslonal parameters 
and fluid property terms. The foundations to the theory are 
presented by Roll and Wislicenus [32] and in a discussion by Aco,feta 
and Parkin [33] of that paper ^ Billet [17] and Fricks [18] 
successfully correlated temperature depression data using this 
method 


18 


Consider a vapor region in an otherwise liquid flow If this 
cavxty is assumed to be steady, then it requires a constant mass 
flow rate of vapor to sustain it. The vaporization Is assumed to 
occur near the leading edge of the cavity at the cavity wall, with 
condensation occurring in the cavity closure region. Thus, the heat 
required to produce this mass flow rate of vapor is given by 

1 . j 

q = X (7) 


or, expressing it in terms of the volume flow rate (Q) , 


q = A Q 


( 8 ) 


New, defining a flow coefficient (C^) as 


= - 2 . 


^ V 


( 9 ) 


The heat required for vaporization becomes 


1 ■ fv * 


( 10 ) 


following the convention of convective heat transfer, an average 
heat transfer coefficient (h) can be defined as 


h = 


- V 


9 


( 11 ) 
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where is the characteristic cavity surface area. Substitution 
of Equation (10) Into Equation (11) and rearranging yields 


AT 


= Di . 




( 12 ) 


or, rearranging in terms of nondimens ional parameters 


AT 


Pe X 
= Nu C 


(13) 


where 


and 


Nu 


- — = area coefficient. 


V D 

CO 


h*D 


Peclet number. 


= Nusselt number. 


and the fluid property terms are evaluated at the bulk temperature 
of the liquid Division of AT In Equation (13) by the fluid 
property terms (p/p ) (A/Cp ) produces the Jacob number on the left 
hand side of this equation. 

The Peolet number and fluid property terms on the fight hand 
side of Equation (13) are known from the free 'stream conditions 
(T_ and V^) and characteristic model dimension (D) - How'ever, C^, 


Cq, and Nu are characteristic of the cavity flow and are to be 
determined empirically. 
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Ihe flow coefficient; (C^) is determined from ventiiaced cavx;^/ 
tests. Ventilated cavities have been used cc model natural cavities 
since the work of Reichardt [34] showed that the drag and geometry 
of ventilated and natural cavities were the same rcr a given cavitation 
number and flow geometry. Based on this result, it is assumed that 
the flow’ rate of vapor required to sustain a natural deveicped ca/ity 
can be estimated by measuring the volume flow rate ot gas needed to 
sustain the corresponding ventilated cavity. This flew rate ot gas 
is expressed as a nondimensional flow coeificient. 

Billst and Weir [26-2'J correlate flow ccefflcienc data with the 


: elation. 


S ^ ill 


U4) 


where 

V *D 

CC 

Re Reynolds number, 

"l 

V 

ot 

^ = Proude number , 

‘^gD 

and 

^ = dimensionless cavity length 

This IS similar tc the ccrrelatlon made hy Reicha>dt [34], using 
the Proude riumfaer and maximum ca^^ity dlametef 

L .. . ; , 

The area coefficient was determined from phocogfmpfis or the 
developed cavity conccur Data by Wald [25], and Billet, Hell, and 
Weir [24] show that the cavity geometry depends only on 
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dimensionless cavity length. Therefore, the area coefficient data 
can be correlated by 




( 15 ) 


Finally, the Nusselt number is determined from Equation (13) 

using the correlations for C. and C_ and measured values of AT- The 

A 

Nusselt number (Nu) is assumed to be of the form 


e f h 

Nu = Re Pr Fr 




(16) 


where 


Pr 


a. 


= Prandtl number. 


Equation (16) employs the standard heat transfer parameters, 

Reynolds number and Prandtl number, plus the important cavity flew 
parameters, Froude number and dimensionless cavity length- 

Once the C^, C^, and Nu expressions have been obtained, 
temperature depressions may be calculated from Equation (13) Or, 
by combining Equations (13) to (16) , AT can be correlated directly as 

r::j 


AT = C, • Pe • 
4 


Ll-^ „ k -r, m V A 
— Re Ft Pr — 7 T" 

^dJ Pj^ Cp 


( 17 ) 


The entrainment theory is certainly not lacking in assumptions 
or empirical constants. However, it is derived from basic concepts 
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and gives a sound physical basis to the problem. It can be used 
successfully to correlate temperature depression data as shown in 
Chapter IV. 


33 B-Factor Theory 

I 

The B~f actor theory is a quasi-static model for the thermo- 
dynamic aspects of a developed cavity. It is assumed that the flow 
is steady, irrotational, and frictionless, that the vaporization 
process is one-dimensional and adiabatic, and that the fluid is 
continuously in stable thermodynamic equlllbriumr Acosta and 


Hollander [4] state that the B-factor theory is artalbgb'us to the 
problem of an insulated cylinder with a tightly fitted piston. The 
cylinder is initially filled with a saturated liquid. Then, as the 
piston is withdrawn, a certain volume of vapor is formed- The B- 
f actor Is defined as the^ ratio of the volume of vapor to the volume 
of liquid involved in the vaporization i^rocess, 

! f ' ^ 

The classic derivation of the B-factor theory is usually 


attributed to Stahl and Stepahoff [ll; however , Jacobs [3] and 
Fisher [2j have also obtained similar results. Assuming constant 
fluid properties, the change In specific enthalpy of the liquid 
phase Is given by 


Ah^ = AT . 


If thermodynamic equilibrium is reached after vaporization, then the 
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energy balance per unit mass of liquid is 


A 


( 19 ) 


But, from the definition of the B-factor 


V- p 


Py 


( 20 ) 


or substituting Equation (18) into Equation (20) 


Pl 

B = at 

Pv ^ 


( 21 ) 


Equation (21) relates the B-factor to the cavity; temperature 
depression t Applying the Clausius-Glapeyron relation for AT and 

Pv 

assuming that — << 1 gives an alternate expression for the B-faccoi 
Pt 


B = 


Pv^ 


C„ T AH„ 
00 V 


( 22 ) 


where AH„ = AP„/p, is the cavity pressure depression expressed as a 

V V Li 

change in head 

! : 

Gelder, Rugger i, and Moore [9]| refined the ealculation oi cbe 
B-factor to account for changes in fluid properties with temperature 
and changes in enthalpy of the vapor r. Again, considering the 
quasi-static model of Acosta and Hollander [4], the piston is 
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withdrawn in small Incremental steps such that the system remains 
in thermodynamic equilibrium« The energy balance (19) per unit 
mass of liquid at some step n is given by 




X = 


n-1 

1 - - E My 

n i=l i 


n-1 

Ahj^ + E My Ahy 

n 1-1 n n 


(23) 


where the effect of the change in enthalpy of the vapor has been 
added, and the summation accounts for the vapor created in previous 
steps- The volumes of liquid and vapor are given by 


n V. 


= E 


V 


i=l 


n 


(24) 


and 




n 

1 - E My 

i=l i 
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io that the B-factor becomes 


B = 


V 


n 

E M^. 

i=l i 


n 


1 - E 




"1 -'Pij ' 
I ! V| 


i=l '1^ 


n 


(26) 


where n is the number of incremental steps and the My are 

1 

calculated from Equation (23), Note that if n = 1, Equation (26) 
reduces to Equation (20) . 

Due to the lengthy calculations required for the Gelder, 
Ruggeri, and Moore [9] method, Hord and Voth [12] devised a 
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simpler method' based on entropy considerations co Account for fluid 
property variations in the calculation of the B-factor. Since, from 
the assumptions of the quasi-static model, the vaporization process 
is isentropic, then, as shown in Figure 36, the initial (i) and 
final (f) entropies can be equated as 


”i. “l. - V 

11 f f f f 


( 27 ) 


and conservation of mass gives that 


\ + V 

iff 


<'28) 


Combining Equations (27) and (28) with the definition of che 
B-f actor yields 


B = 




Li 


IPyi 

f f ^ 


(29/ 


The B-faccor may be calculated from Equations (21), (26), or 
(29), once the temperature depression is known However, none of 
cbese equation-s-a^'r.ount for the dynamics of the cavity riow. 
Gelder, Ruggeti, and Moore [9] show that the B-faetcr should be 
proportional to the maximum cavity diameter (D^^) and inversely 
proportional to the thickness of the liquid thermal boundary layer 
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(5^) Since, in general 


and 


= f^(L/D) ^ 


= ^2 


V *D 

00 


(30) 


(3J ) 


the B-factor is correlated by the relation 

n. 


B = B 


R 


■(a,) 1 

i 

r 

i l/d 

"2 

"a 

fD_| ^ 

V 1 
1 

^ • 


ia/D)^j 

1 

Id 1 

V 1 

CC 1 

R- 


(32) 


where the subscript R indicates a reference value Hord [14-161 has 
inttcduced additional scaling parameters such as surface tension, 
kinematic viscosity, and two phase sonic velocity 

Thus the procedure used in applying the B-factcr theory is tc 
calculate the B-factor for a measured AT from an equation such as 
Equation (29), and correlate the results in a form such as 
Equation (32) where one data point is used for the reference values 
Then, temperature depressions may be estimated by working backwards 
through the procedure. 

The B-factor theory and the entrainment theory ate both 
seml-empirlcal in nature. The entrainment^ theory, based on a 
dynamic model, is more closely related to the actual flow field 
than the quasi-static B-factor model. However, by comparing 
Equations (13) and (21) it can be seen that the B-factor theory 
and the entrainment theory correspond if 
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B 


C. Nu 
A 


(33) 


Thus, employment of the entrattnnient theory enables the B-factor to be 
expressed in terms of the four dimensionless parameters which describe 
the geometric, entrainment, and heat transfer characteristics of the 
cavity 

3 = 4 Two-Phase Laminar Boundary Layer 

In many cavity flows the length of a developed cavity is much 
greater than its thickness, and the flow resembles that of a boundary 
layer or wake. The solutions for such flows are well known in 
this section, methods developed for the solution of one phase 
boundary layer problems will be applied to a two-phase flow 

Consider a heated, two dimensional, infinite flat place as 
shown in Figure 37 with an incoming flow velocity and temperature 
T^ The plate heats the liquid so that a vapor layer is formed on 
the place. The saturation temperature of the liquid is assumed to 
be the same as the plate temperature T, so that the vapor is 
isothermal The flow is assumed laminar with constant fluid 
properties, and the boundary layer approximations are assumed to be 
valxd in the vapor and liquid region near the plate- This model 
ccrresponds to an infinite cavity with zero cavitation number 
However, heat is convected away from the "cavity" in this problem 
(T^ - T^ < 0), whereas in the natural cavity case, heat is 
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convecced to the cavity > 0) , but the heat transfer 

mechanisms should be similar* 

Refer to the coordinate system in Figure 37 The conservation 
equations in boundary layer approximation for the liquid are 


n ‘ • *. 8u 3v _ 

Continuity: — H = 0 , 

3x 3y 


(34) 


Momentum: p 
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Energy: 
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and for the vapor are 


Continuity : 


Momentum: p 


V 


and 


3u 
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3u 
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3x 
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(38) 


Energy: T = T = constant , 


(39. 


where the tilda indicates dimensional variables. The boundary 
conditions for the problem are 


u(x,Q) = v(x,0) = 0 , 


(40) 
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T(i,6) = , 


(41) 


and 


u(x,y) = as y -> “ 


T(x,y) = T as y 


(42) 


(43) 


Three conditions are also required at the liquid-vapor 
Interface (y = 6) ^ Velocity continuity states that 


“liq ■ Sap “ ? ‘ * 


(44) 


and stress continuity requires that 


3u 




= U, 


LIQ 


8u 


3y 


at y - 6 


VAP 
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A third interface condition can be derived from mass continuity f'x 
a control volume as shown In Figure; 38, 


= M „ + A „ , 

AB BD AC CD ’ 


V.46) 


since the flow is steady ^ These mass flow rates can be expressed 
in terms of velocities as 


"ab - h " T '** • 

dx 


( 4 ?) 
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= Pv ^ ’ 
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and 


\c = Pl 


Py 


V + — dy 

3y 


~ 3u - 
u 4 dx 
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dx 


dx 


(49) 
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Thereforej neglecting terms of order (dx) , the mass ccncinuicy 
relationship is 
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The problem can be simplified by introducing the dimensicnies: 
variables 


(52) 


and 


X = 


X 

D 




T - T 


0 = 


- T 
C “ 


(53) 
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where is the characteristic; velocity scale and D is a chata:ter- 
istlc length scale Also, to satisfy the continuity relationf 
Equations (34) and (37), stream functions cj) and ){/ can be 
introduced such that 
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= 


dy ’ 


_ li 

9x 


in the liquid , 


u = 


and V = - in the vapor 

9y 9x 


Therefore, the governing equations become 
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and AT = T^ - T^. The boundary conditions become 


i|j(x,0) = 0 , 


^ (K,0) - ^ 


9x 


9y 


(x.O) =0 , 
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Expressing the total heat transferred (qrj,Q,j.) in terms of a local heat 
transfer coefficient (h) and using known expressions for and 

^COND’ energy balance is 


• 9 T 

h*AT = A - k. — 


9y 


y=6 


(71) 


From the right hand side of Equation (51), 


\ = P 


V 


~ d6 

u V 

dx 


which, expressed in terms of the stream function in dimensionless 
form, is 


\ = P„ u 


V OO 


[9y dx 9x 


(72) 


Substituting into Equation (71) and expressing the results in 
dimensionless form, the interface energy balance is 


hD ^ Pe A 

kL ■ "" ■ Pl S. 


l_9y dx 9xJ 


li 

9y 




(73) 


where Nu^ is the local Nusselt number » This completes the closure 
to the problem 

Determination of an exact ^solution to this problem is a 
formidable task- However, an approximate solution can be obtained 
in the limit as p^/p^ approaches zero. For most fluids at a 
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temperature well below the critical temperature, the ratio of vapor 
density to liquid density is small. For example, in water at 38^C, 
the saturation density ratio is 0.461 10 This limit 

15 equivalent to taking the limit as the vapor Reynolds number (R^) 
goes to 2 erc and the liquid Reynolds number (Rj^) goes to infinityn 
for example, tcnsider water at 38'C with = 30^5 m/sec^ , and 
D = 0 655 om Then R^ - 2.78 x 10^ and R^ = 89.1= Thus, it can be 
seen that this approximation is fairly goodc 

The apprcxlrnate solution for the stream function in the vapor 
phase will now be determined. Equation (57) can be rearranged in 
the iorm 


R R ll- 

V 9y 9x ^ 9x 9y ^^2 


Ml 

8y 


= 0 


(74) 


In the limit as approaches zero. Equation (74) becomes 


9y" 


(75) 


>£ upon integrating three times. 


ip := c^y^ + C^y + , 


( 76 ) 


where C,, and C.. are functions of x only. Then, applying the 
boundary condlticns Equations (62) and (63), Equation (76) is 
reduced to 
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i}j = C^(x) 


(77) 


where Cj^(x) must be determined from a matching condition at the 
interface 

Similarly, the approximate solution for the stream function in 
the liquid phase can be determined By taking the curl of Equation 
(58), It can be con’werted to the vorticlty transport equation, which 
in boundary layer approximation is 


[8y 3x 3x 3y 


= 

9y" 


(78) 


In the limit as approaches infinity. 


j 34 3 34 3 

1,^ 3 ;^ ■ 3l ^ 


lli = 

3y2 


(79) 


A iirst integral of Equation (79) shows that the vorticity is 
constant alcng a streamline- Since the flow upstream of the plate is 
irrotaticnil , then the vorticity is zero everywhere in the liquid and 


9y' 


(80) 


Equation (80) can be integrated twice to yield 


4 = K^y + 


(81) 
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As a result of Equation (64) Kj^ = 1- The function K 2 must be 
determined from the matching conditions at the interface » Therefore 

I . . . ' 

(J) = y + K 2 <x) . (82) 

There are two unknowns in Equations (77) and (82) and three matching 
condxcrons. The problem is over-constrained since the order of , 
Equation (78) was reduced when the limit as approaches infinity 
was taken. Equation (68) leads to a trivial solution for ip 
(C^Cx) = 0) and therefore this matching condition will be 
discarded Physically, this decision neglects the condition of 
stress continuity, which determines the jump in slope of the 
velocity profile at the interface (see Figure 37), in favor of the 
more important requirement of velocity continuity across the 
Interface. The remaining interface conditions are 


3y 3y 


(83) 


and 



at the Interface (y = n^) . The final form for the vapor stream 
function is determined by substituting Equations (77) and (82) into 
Equation (83) and evaluating them at y = Therefore 
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(85) 


Siniiiarj.y, the final form of the liquid stream function can be 
determined by substituting ijj and (j) into Equation (84) and taking 
the limit as P^/P-^ approaches zero« The resulting expression for 
(x) is 


dx dx 


( 86 ) 


or , 


^2 ■ *^6 S 


(87) 


The constant can be arbitrarily taken as zero, so finally 


4) = y - hg 


( 88 ) 


and the =pprcximate expressions for the stream functions have been 
determined- 

There is one assumption in the above derivation that has not 
been explicitly stated- The interface was assumed to be a streamline 
for the flow in the liquid phase when the limit as Py/Pj^ approaches 
zero of Equation (84) was taken. However, there is a velocity 
component normal to the interface, otherwise the vapor layer could 
net grow along the plate. Therefore the interface is not really a 
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streamline, but if the density ratio is small the 'leakage' across 
the interface will be small. 

The energy equation in the liquid was given by Equation (43) as 


3(j) 39 3(|) 36 _ 1 

3y 3x 3x 3y Pe 


3y 


2 AT 


3‘^e 


|i!l 

^3y^i 


(89) 


But $ = y - and taking the limit as Pe approaches infinity. 


30 ^ 30 ^ 

3x dx 3y 


(90) 


The sclution to Equation (90) as obtained by separation of variables 

IS 


0 = C e 


K(ng-y) 


(91) 


where C and < are constants The boundary condition Equation (66) is 
trivially sarisfied and Equation (65) gives C = 1 Thus 


0 = e 


K(n^-y) 


(92) 


where < is as yet undetermined. An expression for the local Nusselt 
number can now be determined by substituting the approximate 
sc/iutions fcr \p and 6 into Equation 03), The result is 


Nu = 


Pe A 


1 


X AT p, G_ 1 2 3x j 


+ < 


(93) 
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Here < represents the contribution to the overall local Nusselt 
number due to the effect of the liquid conducting heat away from 
the vapor layer If it is assumed that the liquid conduction effect 
is small, then K can be taken as zero. Solving for AT, 


AT 


Pe Py A [l 

Nu^ Pl C 12 J 


(94) 


The term 




is the nondimenslonal volume flow rate per unit area 


across the interface. This corres 


- , Q 

ponds! directly to the term ^ in 


the equation for tljie entrainment theory , Thus by comparing Equations 
(13) and (94) , one can see chat this approximate solution to the 
ccnservarion equations produces an expressxon similar to that 
employed in the entrainment theory. 


3 . 5 Nusselt Number Estimate from a Turbulent Boundary Layer Model ‘ 
In the entrainment theory development, an empirical relation 
was assumed for the Nusselt number using the relevant dimensionless 
parameters- in this section a theoretical estimation of the Nusselt 
number is made by applying a caviry model proposed by Brennen [30]. 
This cavity model has been success fully used by Billet and Weir 
[26-2^] to accDunt for gas diffusion effects on ventilated cavities 
Following the procedure outlined by Brennen [30], the flow over 
a cavity wall is represented by flow over a flat plate. As shown 
In Figure 39, the x-axis is placed in the direction of the flow, 
the y-axis is perpendicular to the flow, and the turbulent boundary 
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layer has an approximate origin near the cavity formation points 
Initxally, the velocity and temperature profiles in the liquid are 
different, however, it is assumed that the profiles will be similar 
at the end of the cavity provxded thdt the turbulent Prandtf number 
xs nearly unity. This change in the temperature profile is caused 
by the flux of heat into the cavity. 

The difference between the temperature profile at any position 
X and the xnxt rally uniform temperature profile represents a loss 
of heat from the boundary layer which can be written as 


Aq = 


Pf Cp 


L *' 


(T^ - T)S dy . 


(95) 


Sxnce the velocxry and temperature profiles are assumed similar, 
they can be expressed as defect profiles in the form 


U - u T - T 

CC 00 

U - U f ^T^ f (y) . (96) 

oo C “ C 


where and T^ are the velocity and temperature at the cavity 

I'"' 

surface and and T^ are the velocity and temperature outside of 
the boundary layer- Substituting Equation (96) into Equation (95) 
and evaluatxng at x = L gives the total heat flux into the cavity as 




Pl 



U AT 


00 1 


r' 

U - U ' 


f(y) 

1 - 

00 C 

u 

f(y) 

0 

— 

00 

« t 

. — 


dy 


(97) 
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where is a characteristic cavity size parameter for an 
axisymmecric cavity and AT = 

Since the cavity surface is assumed straight, the turbulent 
boundary layer has the same boundary conditions as half the wake flow 
behind a flat plate parallel to the free stream, provided that the 
ratio of momentum thickness (0^) to cavity length (L) is small. Then 
from Townsend [35], f(y) is a Gaussian function given as 


f(y) = exp 





IT 


2 


) 


^t' 


(x-x ) 
o 


(98) 


where the momentum thickness is assumed constant along the free 
streamline In Equation (98), x^ is the virtual origin of the 
turbulent boundary layer, R^’ is a constant representative of the 
large eddy structure, is the drag coef flcie^nt, and D is a 
characteristic dimension of the body^ 

The value at the drag coefficient can be shown to be 
pr opc r tional to the momentum thickness at cavity formation (6jjj) 
ana inversely proportional to the body size so that 


C_ = constant • 0 /D 
D m 


(99) 


The length of the turbulent boundary layer to the end of the cavity 


(x-x^) can be approximated by the cavity length, which together with 
Equation (99) gives 
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t(y/ = exp 

Mow, cobst iccting EquaLlcn ( 100 ) ,inc.o Equatico ( 97 ) and neglecting 
re.*m5 cl a^der yxelds 

X — . 

• [ I 

q - -cnsLant • D p C U dT = exp - 7^— (dy (lOlj 

^ ’ L 'Ll 

^f.tegiaticn or chis equa'ion gives che resa_>.L '.hat 

q - --.ns rant • hi v'6^*L , (102; 

1 . 

whc : e h“ m.mencurri rhvckness (0 > and rhe ;,har*: ceristic diamerer of 

m 

't- ;av:i.'> axe , oe evaluated for cbe parricular mcdel geomeriyj 

•:.c d - r.e cursiae the boundary layer is assumed equal to "he 

tree si'cair, •-...'icy can be noted that oegietting the higher 

eras' 'r'fr.s it rriw- case -s essentially chs same apprcximaticn as 
p'.du-ea in 'hs -dminar b...ndary layer de taipment in Section 3 A- 
The ^^-^..^0. -r dsSumed to na/e a constant -elcricy in the y-dlrectlcn, 
bu". f.a-e s': in.'i van at ten in teraperat.. ' e However, the 

'Sipera' ..r'^ p.v-iilss a*= diii = :en". t i the r»; rases 

f.t ctfc axisymmer . : c igive bodies considered in this study, 

"ha ra life : IS : 1C oavity diameter (Dj^; can be taken to be 
pr up . t . .nax r: the model diameter (D) and the mcraentum thickness 
e^iimateci from 'he txS' plate laminar boundary layer expression 
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0 = 
m 


0 664 2 


fv ‘xl 


0.5 * 


(103) 




where the distance x can be approximated by the model diameter. This 
gi -es an expression for the heat flux as 


q = constant C C AT (L/D)° 0*75 0,25 _ 

L OD L 


Expressing this in terms cf the Nusselt number. 


Nu = 




ata^X 

W L 


= constant • 




L Pl “ L 


or, in terms of dimensicniess quantities 


,, CCDS rant ,, .,^,0,5 
No =■ r (L/D) 
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CC 1 

L 




(105) 


As determxned by Billet, Hall, and Weir [2‘4], the area coefficient 
the Ee-c-caiiber ogive can be expressed as 


C = constant ' (L/D) 
A 


1 19 


(106) 


and f,c the quarter-caliber ogive, 


1.18 


C = constant ' (L/D) 
A 


(107) 
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Since these two expressions are almost identlcalj the theoretical 


estimace of the Nusselt number for either model is 

Nu = constant • (L/D) ^ Pr v (108) 


This expression will be compared to the empirical determination of 


the Nusselr number in Chapter IV 


CEAPTER IV 


EXPERIMENTAL DATA CORRELATIONS 

4 i I nxcoduction 

Correlation of experimental data is an important aspect of both 
the entrainment and the B-factor theories. In addition, correlation 
of experimental measurements of cavitation numbers, cavity geometry, 
and :avxty liov coefficients help to simplify analysis of these data 
In this haptet, the cur »e-fit method used to correlate the data in 
tb-s study IS presented, and the resulting empirical correlations 
are pcssentea and dis .ussed Data correlations obtained by other 
incest igaccrs are also given fcr comparison 

^ Cufve-Fac Metho d 

In Chapter HI, empirical relations rcr the various 
dimerisionless ca..-ity flow parameters were assumed to be of the 
general form 

F = C - g“ • , , (109) 

I ' ' 1.. 

where F is the quantity tc be ccrrelated, G, H, . are the 
relevant scaling parameters, and C, Q, 6, ate constants- Taking 
the natural lugarithm cf both sides yields 

^ 5 * 
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Ir.f *• inC Q InG ♦ 6 InH • 1 110) 

This equation iinear in terms ot the icgatithm oi the parameters 
The stindard icast-aquare = c^rve-iit method can be applied (see, 
for example, Becket and Hurt (36J> to solve to' the unkncwns C, Q, 

B, . 

The procedure invcKes min:miaing the sum of the squares cf the 
difierence between the experimental and calculated values ci the 
logarithm ot F, which ^an be expressed as 

N 

L ( .fiF - InF)' = /\ , rill) 

where N vS the number ,t data points and the tiida indicates expet i- 
mentai values Subst :iut it,g Equat --n (il0» intc Equation (111) 
y le ids 
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I .aF - r - a *r.u - b Ir.H - )'• a , (Hi:) 
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der; s'^ ..e is t^ken <) 'h •‘-spe. i. eai h unKnowr, , a, B, 

ahd set equ,>; to zero Ih.;= yit,.as a -.et •_! s imu . tane^us equations 
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and 0.318 cm. This corresponds to model to test section diameter 
ratios (D/D^) of 0.167 and 0.083 Thus, the blockage effects fcr 
these data are important. In addition, data were obtained by 
Hord [15] on the quarter-caliber ogive for a D/D^ equal to 0 308 
Details of the analysis of the data obtained by Hord are given in 
Appendix B. The resulting data correlations in the form 

O “ (L/D)" (1141 

are given in Table III for the two model contours and various values 
of D/D^. 

The cavitation number correlations in Table III are plofed in 
Figure 40 for the zero-caliber ogive and in Figure 4i for the quarter- 
caliber ogive The symbols on these figures indicate an average 
of at least 40 data points. Referring to the figures, it can be 
seen ••hat as the model to test section diameter ratio increases, be 
cavitation number increases for a given dimensionless cavity length 
Thus, blockage effects can significantly affect the cavitaticn number 
In fact, the data by Hord [15] is so strongly affected by bio;kage 
that he assumed that the cavitation number was independent ct 
dimensionless cavity length and grand averaged all of his data 
H..wever, the curve-fit of his data in this study does show a weak 
dependence on L/D as shown by the upper curve in Figure 41. 


I 


One method to correct measured cavitation numbers affected by 
blockage is to determine the blockage correction factor (N„) defined 

D 

as 


^ ^ cavitation number (115) 

B unblocked cavitation number 

Figure 42 gives the blockage correction factor as a function of 
dimensionless cavity length and model to test section diameter ratio 
for the two model contours. It can be seen from the figure that N 

o 

increases with increasing L/D and D/D^ The blockage correction 
factor also appears to be nearly independent of nose contour for the 
models considered. Expressions for N may be determined from the 

O 

corresponding cavitation number correlations in Table III. 

4 4 Entrainment Theory Correlations 

As discussed in Chapter III, employment of the entrainment 
theory requires empirical correlations of the cavity area coefficient, 
■ avity flow coefficient, and Nusselt number The correlative 
expressions are given by Equations (15), (14), and (16) respectively 
The correlation of the cavity area coefficient was made from 
the data given by Billet, Holl, and Weir (24) for zero- and quarter- 
caliber ogives with model diameters of 1.27 and 0 635 cm for a range 
cf velocities and dimensionless cavity lengths in water These 
correlations are shown as the solid lines In Figures 9 and 10 
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The correlations of flow coefficients were maae from the data 
given by Billet and Weir [26-27] for zero- and quarter-taliber ogives 
with model diameters of 1,27, 0 635, and 0 318 cm tor a range ot 
velocities and cavity lengths In water These correlations are 
shown as the solid lines In Figures 11-16 

Correlations of Nusselt numbers were made trom maximum 
temperature depression data obtained in this study tor zero- and 
quarter-caliber ogives with model diameters cf 0 635 cm in water and 
model diameters of 0.635 and 0.318 cm in Freon 113, and data obtained 
by Fricks [18] for zero- and quarter-caliber ogive; with model 
diameters of 0,318 cm In Freon 113 The ranges ct velocity, 
temperature, and dimensionless cavity length employed in this study 
are given In Table II, In addition, a correlation for the Nusseit 
number was determined from the temperature depression data taken by 
Hord fl5] on quarter-caliber ogives in nitrogen and hydrogen (see 
Appendix B) 

The constants and exponents for all of the data correlations are 
given in Table IV. These include the correlations the zero- and 
quar t e r -caliber ogive data obtained in th-s study, the correlations 
lot the quarter-caliber ogive data obtained by Horn [15], and the 
theoretical Nusselt number relationship obtained in Chapter 111 
The pertinent equations are referenced in the table. 

Once the correlations of the area coefticient, flow coefficient, 
and Nusselt number have been determined, an expression for the 
maximum temperature depression can be determined from Equation (13) 
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The resulting correlations In the form of Equation (I"’; are also 
given In Table IV and plotted as solid lines on Figures 30 to 32 for 
the zero-caliber ogives and Figures 33 to 35 for the quar ter -caliber 
ogives 

Comparison of the correlative expressions can be made by 
referring to Table IV. The cavity area coetficienr (C^) tor the 
two model contours have almost the same dependence on dimensionless 
cavity length (L/D), and C for the zero-caliber ogive is about twice 

t\ 

that for the quartet -caliber ogive The cavity flow coetficien': 
correlations also have about the same dependence cn ID However, 
the Reynolds number (Re) and the Froude number (Fr) exponents on the 
Cp correlation for the quarter-caliber ogive are greater than fcr the 
2€tc-caliber ogi>/e, reflecting the greater dependence of on 
velocity and model diameter for the more streamlined body 

Comparison of the three empirical Nusselt number (Nu) 
expressions indicates relatively good agreement except for the 
Prardtl number (Pr) exponent The L/D exponents range ticm - 0 ’4 
'o - 1 35, the Re exponents from 1 00 to 1 39, and the Ft exponents 
1 1 ora 0 '’r* to 0 56 All of these exponents have the same sign and 
ate of the same order for like terms However, 'be exponents on Pr 
vary greatly from 0 05 tc 0 84 Thus, the role or the Prandtl number 
in the entrainment theory correlations is unclear 

The theoretical estimate of the Nusselt number derived in 
Chapter 111 is given in Table IV for comparison with the empirical 
correlations It can be seen from the Table that the theoretical 
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exponents on L/D and Pr are slightly greater and the exponents cn Re 
are slightly less than the range of exponents determined from the 
data However, the overall agreement between theory and experiment 
16 quite good. 

Examining the maximum temperature depression (AT ) 

maX 

"orrelat ions, it can be seen chat the exponents all have similar 
signs ter like terms, and they all agree qualitatively except tor the 
Prandti number exponents. Thus, although the data let the 

iero-caliber ogives shows a trend with velocity and model diameter 
which is opposite that of the quarter-caliber ogi e data, the 
AI expressions are consistent when expressed in terms ci 

[TtaX 

dimensionless parameters in that the exponents on iike terms have the 

same sign. In addition, referring t-'' Figures 30 to 35, the empirical 

'T correlations are a goed curve-tit to the experimental data, 
max 

with the standard deviation being 0 25 C tor the aerc-caliber cgive 
and 0 07 C for the quarter-caliber ogive 

4 5 B-Factor Theory Correlations 

Correlations of temperature depression data i:r de.elcped 
ca/i*^ies have been made by several Invest Igat : rs using the B-fa:tot 
theory A summary of the various correlations are given in Table V 
As shewn in the table, three different methods ot calculating the 
B-iacror have been used, and the correlations are expressed in the 


term ot Equation (32) 
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The exponents on the liquid thermal diffuslvlty (a^^) vary 
widely from 3.52 to - 1.00 showing no consistent trends This is 
similar to the problem with the Prandtl number exponents in the 
entrainment theory. The model diameter exponents are also 
inconsistent, but the only two experiments that varied model diameter 
had very different model geometries. The free stream velocity (V^) , 
cavity length (L) , and kinematic viscosity (Vj^) exponents are each 
fairly consistent. Direct comparison of the B-factcr and entrainment 
theory correlations were not made due to the fact that different 
physical reasoning was used to develop the two theories, and 
different curve-fit methods were used to determine the exponents 


CHAPTER V 


SUMMARY 

5 . 1 Theoretical Considerations 

The entrainment theory and the B-factor theory are semi- 
empirical methods to describe thermodynamic phenomena in developed 
cavitation They have been successfully applied to correlate 
measured cavity temperature depression data The entrainment theory > 
however, is a more physical approach derived from a simple energy 
balance that accounts for the cavity geometry, vapor mass flow rate, 
and heat transfer mechanisms. The entrainment theory can also be 
derived from an approximate solution to the equations of motion 
for a two phase, two-dimensional laminar boundary layer A 
theoretical analysis of the cavity heat transfer mechanism based on 
a method proposed by Brennen [30] yields an expression for the 
Nusselt number which agrees in part with empirical results 

5 . 2 Cavitation Numbers and Blockage 

The cavitation number based on cavity pressure is an important 
parameter for describing developed cavities The experimental data 
show that the cavitation number Is solely a function of dimensionless 
cavity length and is independent of velocity, temperature, model 
diameter, and working fluid for a given model contour and model to 
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test section diameter ratio. The cavitation number, when based on 
a local pressure within the cavity. Is nearly constant for the first 
60% of the cavity, then decreases rapidly near the cavity closure 
point due to the Increase in cavity pressure 

Blockage can have a significant effect on measured cavltaticn 
numbers. A good measure of tunnel blockage for axisymmetric flows 
is the ratio of model diameter to test section diameter. As the 
tunnel blockage increases, the cavltaticn number also increases for 
a given value of dimensionless cavity length. For severely blocked 
conditions, the cavitation number can appear to be independent of 
dimensionless cavity length. The blockage ccrrection factor, 
defined as the ratio of the cavitation number with blockage to the 
unblocked cavitation number, is an effective means to correct for 
blockage effects in that it is nearly indepenaent oi nose contour 
Empirical correlations of the cavitation number aod the blockage 
correction factor as a function of dimensionless .a\)i>^y length 
have been obtained. 

5 3 Cavity Area and Flow Coe f ficients 

The cavity area coefficient and cavity tlow coefficient are two 
important dimensionless parameters employed in the entrainment 
theory The area coefficient, a dimensionless representation of 
the cavity surface area. Is determined from photographs cl developed 
cavity profiles. The flow coefficient, a dimensionless 
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repi esentat ion of the volume ftow rate ot vapor through the cavity, 
is determined by measuring the volume tlow rate ct gas needed to 
sustain a ventilated caxity. 

The experimentally determined area coefticient is solely a 
function of dimensionless cavity length For different model contours, 
the magnitude of the area coefficient changes, but its dependence on 
dimensionless cavity length remains nearly the same The 
experimentally determined flow coefficient increases with increasing 
velocity, model diameter, and dimensionless cavity length For 
different model contours, the dependence of the flow coefficient on 
dimensionless cavity length remains nearly the same, but its 
dependence on velocity and model diameter is greater tor the mote 
streamlined quarter-caliber ogive than tor the ^e.' 0 -callber ogive 
Empirical correlations of the area and flew ccetficcents with the 
pertinent tlow parameters have been determined 

5 ^ Measured Temperature Depressions 

Measurements of the difference between the tree st'eam 
temperature and the cavity temperature were made icr a variety of 
flow conditions. These data, along with data obtained by other 
.nvestigators, provide a wealth of information cm erning thermo- 
dynamic effects on developed cavitation 

The axial variation oi temperature depressions within the 
..avity was found to be linear, with the maximum temperature 
depression occurring near the leading edge ot the cavity. The 
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nidXiiT.um temperature depression increases with increasing 
dimensionless caviry length and tree stream temperature However, 
the maximum temperature depression increases with in.reasing tree 
stream velocity and mcdea diameter t^r the quartet- aliber cgives, 
but has the opposite trend for the zero-caliber ogives Ali .iT 
data show consistent trends when analyzed in teims ot Reynolds 
number, Pranatl number, Pec let number, Froude number, and 
dimensionless cavity length :n the ctntex» oi the entrainment fheoiy 
because the signs ol like terms ate the Same The maximum 
tempera'uie dep;ession data have been used to aetermite Nusselt 
number relations fer aeveloped cavitiea, whi.n, when employed with 
the entrainment theory, can be urtd i.u estimate ta itv temperature 
depressions . 

5 5 Re-v.m met.dat icns Icr further St udy 

Three basic pr'^blems .on.erning r he rmodvnamic streets on 
develvjptd cavitaticn need lurthtr study first, the titter ot mode! 
coniigurati n on cempetature deprtsti-ns has n..i bett. adequar.e*v 
determined Additional expet imtn* la data is needed tc extend the 
applicaCi.n Ol the ent'ainmtnt theory to aaai'i.aol t«_ and thi.ee 
aimensi->naa bodies 

Second, an extension ot present theciier. ^,hcula be aitemptea 
»hat w„uad account ter the eitects ot model gc_mttry on the cavity 
heat transfer mechanism^ In addition, m.rt ex'tnrive and general 
rheors'i.al analyses ct thermodynamic ette-ts should be conducted 



Third, the role of the fluid property terms, especially the 
Prandtl number and/or thermal diffusiviiy, in the thermodynamics 
of developed cavities remains unclear. Further experimentation 
IS needed ever a wider range of temperatures in many different fluids 
to further define the effects of fluid property variations Also, 
attempts should be made tc determine it vapor fluid properties other 
than vapor density have an effect on the ._a.<ity temperature 
depression 
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APPENDIX A 


THERMOCOUPLE FABRICATION AND 
INSTALLATION PROCEDURE 

Cavity temperatures were measured with ..pper-cDnstantan 
thermocouples installed in the cavity temperature models Three 
thermocouples were constructed trom 0 010 ,.m diameter wire sto.K 
with a 0 00025 cm polyurethane coating I;, i.rm the the- rmo.ouple 
bead, the insulation was caretully stripped away trem one end cf the 
wires ana the copper and constantan wires twis'ea together The 
wires -vfcre then welded with a helium atmosphere arc-welder whi'h 
termed the bead. The bead was reiniorced with a small amount ct 
epoxy cement and the leads coated with enamel paint to prevent 

sb 0 r , 

The r he rmocoupies were then string tirjugh the therm,, couple 
ports and out the rear of the test m,dei Ea h str ot *eads was 
Inia 'b'ough 26-gauge teflon -ubing that extended inro the model 
fc: added protection The the rmor Ooples were positioned in the ports 
su:h that the beads ptccruded just abeve the model surface The 
remaining .space in the port was filled with epexy cement to secure 
the bead tn position The tetlon tubes holding the thermocouple 
leads were an bored to the rear ct the mode* using General Electric 
Rr>' compound. The model was then inserted :nto the sting- mount 




itiiiii-rn~"niriTirtTrnTT-. 
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with 'he tubing tunning through one c£ the stings. Finally, the 
tubing with the thermocouple leads inside was anchored inside the 
St ing with epoxy 

The model was then installed in the tunnel, with the leads 
exiting through the runnel wall The three thermc.ouples were each 
connected in series with the tree streair. t herm_ ic uple , and the 
procedure completed 
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APPENDIX B 

ANALYSIS OF THE DATA OBTAINED 
BY J HORD 

Extensive data of measured pressure and temperature depressions 
for developed cavities on quarter-caliber ogives were made by 
Herd ll5j These data were obtained in liquid hvdrcgen and liquid 
nitrogen for model diameters ot 0 533, 0 907, and l 067 cm for a 
variety of flow conditions The data are presented in tabular form 
at the end of his report, making it easy to analyze The tletmo- 
pbysical p»-3perties of liquid hydrogen and nitrogen ate available 
frem me National Bureau of Standards in References (r.lj and (42] 

To determine the cavitation number as a t unction cl 
d mens J rn less cavity length, the cavity pressure measured at the 
•eading edge pressure port, free stream pressure, and iree stream 
•e.-.ity were used The empirical correlation was then determined 
b> :„'ve-t lit .ng the resulting data Three ditie'ent model 
diameters were tested, but only one value ot mcdel t<. test section 
d.ameter ratio was tested since three seated test se^ricns were 
used 

The entrainment theory was also app,ied to the temperature 
depression data Due to the large volume cl data, linear 
ext rapotat ions to the leading edge cf *he ,.a\ity were not made. 


i 
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4.5 


The tempe j'.-e aep-ession measured the ieua:r.g idge therm oapie 
y>ai c.nsidered tc oe the maximum temperatort dtpreesicn The 
Nusselt number ccr.eeatun was then dtte-min&d using 'rtsr- max. mam 
tempeiatuie depression data and tbe t . jvo ctlti.ien' and a'ed 
ccerficient :tia:io;ns .^oi ined in rhls stuav 
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AVAIIABI.i: KXPKRIMKNTAI. DATA OF 1 HF.KM01)YN.\.V!r KKFECTS ON 
IIFVKI.OPKI) CAVITAiK'N 
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Figure 3: Photograph of Test Models for Cavity Pressure 

Measurements 
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Figure 9: Area Coefficient Versus Dimensionless Cavity Length 

for the Zero-Caliber Ogives 
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Figure 11: Flow Coefficient Versus Cavitation Number for the 

0.318 cm Diameter Zero-Caliber Ogive 
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Flow Coefficient Versus Cavitation Number for the 
0.635 cm Diameter Zero-Caliber Ogive 
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Figure 13: Flow Coefficient Versus Cavitation Number for the 

1.27 cm Diameter Zero-Caliber Ogf.ve 
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Figure 15: Flow Coefficient Versus Cavitation Number for the 

0.635 cm Diameter Quarter-Caliber Ogive 
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Figure 17: Cavitation Number Versus Dimensionless Cavity Length 

for the Zero-Caliber Ogives in the 30.5 cm Water 
Tunnel (No Blockage) 
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18: Cavitation Number Versus Dimensionless Cavity Length 

for the Quarter-Caliber Ogives in the 30.5 cm Water 
' Tunnel (No Blockage) 
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Figure 20: Cavitation Number Versus Free Stream Temperature 

for the 0.635 cm Diameter Zero-Caliber Ogive in 
the 3.8 cm Tunnel (D/D.^, = 0.16 7) 
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Figure 23: Cavitation Number Versus Dimensionless Axial 

Position for the 0.318 cm Diameter Zero-Caliber 
Ogive in the 3.8 cm Tunnel 
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QUARTER-CALIBER OGIVE TEMPERATURE: 20-93 deg C 
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FREON 113 ® 36.6 m/sec J ALL TEMPERATURES 
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Figure 25: Cavitation Number Versus Dimensionless Axial 

Position for the 0.318 cm Diameter Quarter-Caliber 
Ogive in the 3.8 cm Tunnel 
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ZERO-CALIBER OGIVE FREON 113 

MODEL DIAMETER: 0.318 cm ©THIS STUDY 1 AVERAGE 

VELOCITY: 36.6 m/sec ° FRICKS (1974)/ OF DATA 
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Figure 27: Typical Temperature Depression Versus Dimensionless 

Axial Position Data for the 0.318 cm Diameter Zero- 
Caliber Ogive in Freon 113 
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QUARTER- CALIBER OGIVE FREON 113 
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Figure 28: Typical Temperature Depression Versus Dimensionless 

Axial Position Data for the 0.318 cm Diameter Quarter- 
Caliber Ogive in Freon 113 
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Maximum Temperature Depression Versus Free Stream 
Temperature for the 0.635 cm Diameter Zero-Caliber 
Ogive in Freon 113 
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Figure 39 


Turbulent Boundary Layer Heat Transfer Model 
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Figure 40: Cavitation Number Versus Dimensionless Cavity Length 

for the Zero-Caliber Ogives in Blocked and Unblocked 
Conditions 
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Figure 41: Cavitation Number Versus Dimensionless Cavity Length 

for the Quarter-Caliber Ogives in Blocked and 
Unblocked Conditions 
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Figure 42: Blockage Correction Factor Versus Dimensionless Cavity 

Length for the Zero- and Quarter-Caliber Ogives 



